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7T° to two-photon decay in lattice QCD 



OS 

o 
o 

(N 

o 

Q 

(N 



in: 



> 

in 

(N 

o 

On 
O 



X 



E. Shintani* 1 ! S. Aoki 2 3 , S. Hashimoto 45 , T. Onogi 1 , N. Yamada 45 (for JLQCD 
Collaboration) 

1 Department of Physics, Osaka University, Toyonaka, Osaka 560-0043 Japan, 

^Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, 

Japan, 

3 'Riken BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA, 
4 KEK Theory Center, Institute of Particle and Nuclear Studies, High Energy Accelerator 
Research Organization (KEK), Tsukuba 305-0801, Japan, 

5 School of High Energy Accelerator Science, The Graduate University for Advanced Studies 
( Sokendai), Tsukuba 305-0801, Japan 

We calculate the neutral pion (n°) to off-shell two photon (7*7*) transition form factor in lat- 
tice QCD. The transition form factor can be extracted from the three-point function of the form 
(axial-vector)-(vector)-(vector) as a function of off-shell two-photon momentum. Since the axial- 
anomaly plays an important role in the n° — > 77 decay process, we employ the overlap fermion, 
which preserves the exact chiral symmetery on the lattice. After extrapolating to the chiral and the 
vanishing photon momentum limit with a fit function based on vector meson dominance (VMD) 
model, we find that the Adler-Bell-Jackiw anomaly is correctly reproduced. 
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1. Introduction 

In this report we present our lattice study of the K° —> yW yM form factor. The photons could 
be either on-shell or off-shell. The goal of this study is two-fold. The first is to reproduce the 
contribution of the Adler-Bell-Jackiw (ABJ) anomaly ^ to the form factor, which is necessary 
condition for futher applications. The second is an application to the estimate of the hadronic light- 
by-light (L-by-L) scattering amplitude, which is important in the theoretical calculation of the 
muon g — 2, y L [§]. Since, in contrast to the vacuum polarization diagram, no experimental 
result is avaiable to estimate the L-by-L diagram, purely theoretical estimate is required, but it is 
challenging ph. Some of the phenomenological estimates is based on the assumption that the L-by- 
L scattering amplitude is dominated by the effective diagram containing yy* — > 71° and 71° — > y*y* 
transition form factors connected by a pion propagator. The results are in the range y L = 
(80 ~ 130) x 10~ n [§ g § 0, §]. Even within this assumption, these are significant model 
dependence. Therefore, model independent calculation of n° — > y*y* transition form factor is a 
valuable step toward the precise test of the Standard Model 

We study the 71° — » y*y* transition form factor in soft Euclidean momentum region. Similar 



attempt is seen in [JlOp , in which the technique of [ |11[ ] for charmonium decay to two photon state is 



applied. Since the axial-anomaly plays a key role, the exact chiral symmetry and flavor SU(2) sym- 
metry are desirable in the lattice calculation. We therefore apply the overlap fermion formulation. 
We use the Nf = 2 dynamical overlap fermion configurations generated on a 16 x 32 lattice at 



j8 = 2.3 [12], which corresponds to the inverse lattice spacing a =1.67 GeV. Topological charge 



is fixed to Q = 0, which induces finite size effect [13]. For chiral extrapolation, we take four values 



of quark mass, m q = 0.015, 0.025, 0.035, 0.05, approximately covering the range of m s /6 to m s /2 
with m s the physical strange quark mass. 

2. 7t° — > YY transition form factor 

In the continuum theory the 71° — > yy transition form factor, f^yyipx^Pi), is defined through a 
matrix element of two electromagnetic (EM) currents V EM between the pion state and the vacuum, 

d 4 xe^ x (n°(q)\V^ M (x)V EM (0)\0) = s, va ppf p P 2 f J[ o rr (p 1 ,p 2 ), (2.1) 

where p\ , p 2 are Minkowski photon four-momenta and q = —p\ — p2 is the pion four-momentum. 
fifiyyipiiPi) can be extracted from the three point function of an axial- vector current and two EM 
currents via 

J d 4 x J jV _,V+Pl> ' ) (0|r{5A 3 (x)V v EM (3;)V' /J EM (0)}|0) 

= - _ * 71 2 ^valiP'lP2f7t Yr(PUP2) + ---- (2-2) 

with the pion decay constant f n = 131 MeV. The ellipsis represents excited state contributions 
which have the same quantum number with a pion. In the Euclidean lattice calculation, only 
negative q 2 region is accessible, and we mainly focus on the small negative q 2 region. Thus, 



in the following preliminary analysis, we assume that the first term of (2.2) gives the dominant 
contribution, and the excited state contributions are ignored. 
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To evaluate the three-point function on the lattice, we use a combination of the conserved axial- 
vector current (A^ v ) and the non-conserved local EM currents (V EMloc ). With the overlap fermion 
formulation the conserved current takes a non-local form A c ™ (x) = £ v z q(y) l a K^ (y, z; x)q(z) where 
K^(y,z;x) is a kernel derived from the flavor non-singlet axial transformation [|l4|]. r" is the gen- 
erator of flavor SU(2) Lie group normalized by trz a T b = 8 ab . The local EM current is defined as 
V^ cEM (x) = Zq(x)Q e y jl q(x) with a quark charge matrix Q e = diag(2/3, — 1 /3). The renormaliza- 
tion factor Z = 1.3842(3) is obtained non-perturbatively fli~5|]. This local current does not satisfy 
the charge conservation and thus unphysical contamination (lattice artifacts) can appear in the re- 
sult as will be discussed later. Since the EM current consists of the flavor triplet and singlet parts, 
the evaluation of the three point function ^c v 3ytocEMy-kx:EM^ re q UU - es [ Q calculate both connected 
and disconnected quark diagrams. In the following analysis, we assume that the disconnected 
contribution is negligible, which is valid in the flavor SU(3) limit. 

We define two three-point functions of type (axial- vector)-(EM vector)-(EM- vector) and (pseudo- 
scalar)-(EM vector)-(EM vector) by 



Gf v V (P 2 ,Q) = £^ ! ' Gv - ! ' P2 ' V (VA cv3 (x)^ ocEM (y)y / l ocEM (0) 

x.y 1 

= tr[T 3 2 e 2 ]( £ e -^-^2ReTr[5,(0,y)5^ A (y,z'|x)^(z / , 3 ;)7v^(^0)7, 1 ]), (2.3) 
G™(P 2 ,Q) = £ e - ! ' ev - I ^(2m,P rot3 (x)C cEM ( 3 ;)\/ / l ocEM (0) N 

x,y ' 

= 2m q tr[T 3 Q2}( £ ^^- i ftv2ReTr[^(0,x)r P (x,y)^(y,y)7v^(j,0)7 M ]), (2.4) 



where () c denotes the connected contraction, while () denotes the statistical average, S q (x,y) de- 
notes the quark propagator and Fp(x,y) = (1 — D 0V /MQ)y=,{x,y). "tr" represents the trace over flavor 
indices and "Tr" represents the trace over color and spinor indices. The Pi 2 and Q are Euclidean 
momenta of two photons and a pion, respectively. On the lattice they are dicretized as 2nn^ / (L^a), 
with the lattice spacing a and the lattice size L x ^ z = 16, L t = 32. The pseudo-scalar density op- 
erator P rot3 (x) = <?(x)t 3 7 5 [(1 - D ov /M )q](x) forms the axial-Ward-Takahashi identity with A^ v3 . 
The reason for considering ( jL4] ) will become clear soon. In the above equation VA CV denotes the 
forward derivative VA cv (x) = L p [Ap v (x + p) — Ap v (x)]. We apply the source method for spatial 
integral over y with P2 = (0,0,0, P 2 t) and v = 1. After making Fourier transformation for t y and x, 
we obtain G^y (P2, Q) and G p ^ (P2, Q) for each pt = 1 ~ 4 and Pit- The remaining momentum Pi 
is determined by the momentum conservation, P\ = —Q — P 2 . 

The flavor non-singlet axial Ward-Takahashi (AWT) identity for the three point functions can 
be expressed as 

Gf v V (P 2 ,Q) = G p ^(P 2 ^) + {(dX° cEM K cEM )c(P2)8\Q) 

+ (Vl ocEM (8y; cEM MP 2 )8 4 (Q) + ■■■ (2.5) 

where the second and third terms are contact terms coming from lattice axial transformation, 
8ayb _ qy 5 y ll T a T b q + qy pi T b X a y 5 {\ -D ov /M )q for (y^cEMybcEM^ These terms do not contr ib ut e 



to fjflyyiPiiPi) as Q = implies Pi = — P2, for which the three point function ( |2.2| ) vanishes due to 
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its Lorentz structure. Since we ignore the disconnected diagram, there may be some contributions 
denoted by ellipsis. Thus, we calculate the both sides independently to estimate the potential size 
of disconnected diagrams. 

The 7T° — > 77 transition form factor is extracted from the three-point function with the follow- 
ing momentum assignments: 

Py = (0,Pi y ,0,0) at n = 3, or Pi = (0,0,P U ,0) at ju = 2. (2.6) 

Then the only nonzero term is the one proportional to Pi a^20 e /ivaj3- Furthermore, these momentum 
assignments suppress a class of lattice artifact due to the violation of Lorentz symmetry and EM 
charge conservation. Since Pi and P2 are orthogonal to each other and have only one nonzero 
component, Lorentz violating terms such as (Pf + P% v )PiaP2p £ iivai} vanish. Therefore we obtain 

p lat (Pi,p 2 ) = Gf v y (p 2 ,e)/(£Pi„p 2/3 ^ va/3 ) = - J?™^ forr(Pi,Pz)- (2-7) 



aj3 



up to contributions from excited states and disconnected diagrams. Since P\i are defined in the 
Euclidean space-time, the on-shell condition is realized only at zero momentum Pi = P 2 = 0. Note 
that we numerically check the consistency of P lat obtained from G^ and G^, v within the sta- 
tistical error. P lat (P ,P 2 ) from (Ob are averaged over the two physically equivalent momentum 
assignments in (|2.6[). 



3. Results 

First we test if the ABJ anomaly [jl]] 



1 



fo rr (0,0) = 7Z z- r , (3.1) 



is reproduced in the chiral and the photon's on-shell limit. The data are extrapolated to the zero 
momentum and zero quark mass limit assuming a fit function. Here we consider the vector meson 
dominance (VMD) model as an ansatz. Since we neglected disconnected diagrams, the EM vector 
current only couples to p meson. The functional form is 

2 

F VMD (Pi,P 2 ;X a ) = - "» 2 X a G v (P h m v )G v (P 2 ,m v ), (3.2) 

with a free parameter X a , which is supposed to correspond to the axial-anomaly, X a = fnfn°yy^-, 0)- 
The vector meson propagator is 

2 

Ttl 

G v (P,m v ) = -tx-T-v (3-3) 
P z + 

where m v denotes a vector meson mass. The vector meson mass is determined independently 
from an exponential fit of a smear-local vector current correlator. Performing chiral extrapolation 
with a linear function, m v = m Q v + c v m\, we obtain m® = 0.798(19) GeV, which reasonably agrees 
with the physical p meson mass 0.775 GeV. Since we encoded the quark mass dependence of 
P VMD (Pi ,P 2 ;X a ) into m 2 n and m v , only one free parameter X a is left for the fit of the lattice data. In 



4 



It to two-photon decay in lattice QCD 



E. Shintani 




(aP 2 f 

Figure 1: Lattice results for F (Pi,Pz) at a fixed (aP\) 2 ~ 0.154 as a function of {aPi) 2 and fit function of 
p VMD (/ , 1) f» 1 ;X a ) (solid curve) and F VMD+ (Pi,f2;^a,C3,C4) (dashed curve) for each quark mass. Thin lines 
denote the VMD (dashed-dotted) and VMD+(solid) in physical point. 




(a) (b) (c) (d) 



Figure 2: Vertecies of 7t — v y*y* decay. The dashed and waved lines show pion and photon propagator, and 
double-line and dotted lines show vector resonance and its excited state propagator respectively. 



Fig. [j], we plot F lat (P\,P2) at (aP\) 2 = 0.154, which is the minimum non-zero value, as a function 
of {01P2) 2 . The fit results at different quark masses are also shown by solid curves, where the fit 
point is taken to be (aQ) 2 = 0.039. From this plot we can see that lattice data and F VMD (Pi ,Pi;X & ) 
are in good agreement at the lowest momentum. We obtain 

X a = 0.0260(6), (3.4) 

which is consistent with the expectation X a = \/(4ti 2 ) = 0.02533 despite of various approxima- 
tions. Probably the main reason is the exact chiral symmetry, because the existence of the conserved 
axial-vector current significantly reduces lattice artifact. 

Going beyond the minimum value of the momenta, the fit curves deviate from the data point. 
To accommodate this deviation, we modify the fit form by incorporating an excited vector meson 
state as 

2 

F WMD+ (P u P 2 ;X a ,c 3A ) = - ™ K 2 X a c 3 G v (P u m v )G v (P 2 ,m v ) 
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+ C4 CB (G v {Pi,m v ) + G v (P 2 ,m v )) + I — C4 



(3.5) 



with mass dependent couplings, C3 = 2 — c° A + m| and C4 = + c™ m|. In this fit form, the 
propagator of the excited vector meson is approximated by a constant G v >(P,m v /) ~ 1 (see Fig. 
We have four free parameters. This parametrization | fl6| ] satisfies the anomaly relation ([Q|). We 
impose a constraint lim m 2^ () (c3 + C4) = 2, so that the transition form factor in the high Q 2 limit 
scales as l/Q 2 as suggested by perturbative QCD [17]. We note that a phenomenological estimate 
leads to (c 3 + c 4 )/2|p he nom = 1.06(13) [JT^] . Setting the fit range to 0.039 < (aP 2 ) 2 < 0.347, we 
obtain the dashed curves in Fig. [[]. Compared with the naive VMD fit (3^2), the agreement of the 



fit with the data in momentum range of (a/^) 2 > 0.039 is improved especially at larger quark mass 
region. Fit parameters obtained using ( |3.5[ ) are C4 = 1.20(21), c"' = —8.4(8.6), c™ = —1.3(5.8) and 
X a = 0.0243(29). Again X a is consistent with the expectation. Furthermore is consistent with 
unity, which implies that the VMD approximation is reasonable and the coupling with the excited 
vector meson state gives a sub-leading contribution. In Fig. [I] we also show the comparison with 
VMD and VMD+ curves in the physical pion and p meson mass. The difference between the two 
lines are almost negligible. 

The construction of reasonable fit function for f^o^y over higher momentum region is impor- 
tant when estimating the hadronic L-by-L diagram, because it is necessary to integrate fjfyty* (P\ , Pz) 
from zero to infinity in both P 2 and P 2 . This study provides suggestive information that the VMD 
model describes the data of fjfi^t '^2) reasonably well, and inclusion of its excited state signif- 
icantly extends the region to larger momentum. 



4. Summary 

We have studied the 71° — ► 77 transition form factor by calculating the three point function, 
(AW), on the Nf = 2 dynamical overlap fermion configurations. The data of /^o^y. is fitted to 
the VMD motivated functions, and in the chiral and the on-shell photon limits the ABJ anomaly, 
/^oyy(0,0) = 1/(4 7i 2 f K ), is reproduced. It is also found that the fit function based on the VMD 
model well describes the behavior of mass dependence at the lowest momentum in which the 
higher resonance contribution is small. Although the present work is feasible study and there are 
many things which have to be understood, e.g. the contribution from higher excited state of pion or 
vector mesons, disconnected diagram effects, finite volume and fixed topology effects, comparison 
with a prediction from perturbative QCD, our result encourages us to apply to the non-perturbative 
estimate of the hadronic L-by-L scattering amplitude. 
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